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Making use of an Ehrenfest-Maxwell relation we show that in Abelian projected 
SU(2), in the maximal Abelian gauge, the dynamical electric charge density gen¬ 
erated by the coset fields, gauge fixing and ghosts shows antiscreening as in the 
case of the non-Abelian charge. 


1 Introduction 


Lattice studies based on Abelian projection have had considerable success iden¬ 
tifying the dynamical variables relevant to the physics of quark confinement. 
There is no definitive way_as yet of choosing the optimum variables, but in the 
maximal Abelian gaugelllla the U(l) fields remaining after Abelian projection 
produce a heavy quark potential that continues to rise linearlyB. Further the 
string tension is almost, but not-exactly, equal to the full SU(2) quantity; 92% 
in a recent study at (3 = 2.5115a. 

This suggests that we may be close to identifying an underlying principle 
governing confinement. All elements of a dual superconducting vacuum appear 
to be presentOu; in the maximal Abeliaingauge magnetic monopoles reproduce 
nearly all of the U(l) string tension Era. The spontaneous breaking to the 
U(l) gauge symmetry Js.^ signalled by the non-zero vacuum expectation value 
of monopole operator M. The profile of the electric field and the persistent 
magnetic monopole currents in the vortex between quark and antiquark are well 
described by an effective theory, the Ginzburg-Landau, or equivalently a Higgs 
theory-giving a London penetration depth and Ginzburg-Landau coherence 
length Enid. 

Central to finding the effective theory is the definition of the field strength 
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operator in the Abelian projected theory, entering not only in the vortex pro¬ 
files but also in the formula for the monopole operator. All definitions should 
be equivalent in the continuum limit, but use of the appropriate lattice expres¬ 
sion should lead to a minimization of discretization errors. 

In Ref. iHl we exploit lattice symmetries to derive such an operator that 
satisfies Ehrenfest relations; Maxwell’s equations for ensemble averages irre¬ 
spective of lattice artifacts. 

The charged coset fields are normally discarded in Abelian projection, as 
are the ghost fields arising from the gauge fixing procedure. Since the remain¬ 
der of the SU(2) infrared physics must arise from these, an understanding of 
their role is central to completing the picture of full SU(2) confinement. In 
the maximal Abelian gauge a localised cloud of like polarity charge is induced 
in the vacuum in the vicinity of a source, producing an effect reminiscent of 
the antiscreening of charge in QCD. In other gauges studied, the analogous 
current is weaker, and acts to screen the source EJ. (This a tentative result, 
however, without the benefit of the refined definition of flux.) 


2 Maxwell Ehrenfest relation 

We first introduce and review the method due to Zach et al. Ill in pure U(l) 
theories. Consider a shift in a U(l) link angle, 9 v (x o) —> 9 v (x o) + 0 s (xo), in 
the partition function containing a Wilson loop source term 

Z W ({0 S }) = J[d(6„+0 S )} sin 9 W e ^E cos V. 


Since the Haar measure is invariant under this shift, Zw is constant in these 
variables. Absorbing the shift into the integration variable and taking the 
derivative gives 


d 


dd s (x 0 ) 


Z w = 


J[dO ] (cos 9w ~ sin6>w /3A M sin 9^) 


cos 6 


= 0 . 


This can be cast into the form 




(sin 9w sin 9^) 

(cos 9w) 


— ed = T 

C/U X,Xw Ul/ ' 


We use the term Ehrenfest-Maxwcll relation because it is the expectation value 
of what is normally a classical extremum of the path integral — an Euler- 
Lagrange equation. If we define flux using sin 6^,, instead of 9 for example, 
then we get a precise lattice definition of current. 


2 




Before addressing the full problem we first generalize from U(l) to SU(2) 
without the complication of gauge fixing, with a shift U fJl (xo) —> U fJl (xo)U s (xo) 

Zw({U s }) = j[d{UU s )} W 3 (U) W 3 = i Tr(UWUUia 3 ). 

The size of the source is irrelevant so we choose it to be the simplest case, i.e. 
a plaquette. We choose the shift to be in the 3 direction 

de to) Zw = °’ = (* “ > 

giving the Ehrenfest relation 

= <W; w = \Tr{uWuu). 

The notation (S)^ denotes an e derivative 0. 

For (3 = 2.5, 0(W 3 (S) M ) = 0.0815(2), and (W) = 0.0818(1), and the 
difference = 0.0003(3); i.e. zero within statistical errors. 

To cast this into the form of Maxwell’s equation we decompose the link 
into diagonal D M and off-diagonal parts 

Up(x) = D^(x) + O^x). 

We then group terms involving the diagonal part into divE and group all terms 
having at least one factor of the off-diagonal link into the current. 

[P(( s )n)w] u=D = -(A nF^w, (• • -) w = ^ ^ \ 
giving the final form of the Ehrenfest relation 

/ a /-I \ / jdyn.\ i t static, c ^ rstatic 

\^-n^nv)w — Wk /wT-h, j o x , xw = -J v 

This then tells us how to choose a lattice definition of field strength that 
satisfies an Ehrenfest relation: 

F, JV = ^-^Tr(D ] D^DDia 3 ) llu . 

The effect of gauge fixing gives 

Zw({U s }) = J [d(UU s )] W 3 {U) A fp <5[F] e^ u \ 
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where we have introduced 


1 = A FP J ndgMYlsmu**™-,*)}, 
3,y 


and integrated out the g variables in the standard way. So A fp = det M 
where 


Mix]jy 


dF?(x) 


dgj{y ) 


3=0 


In this case Zw is not invariant. The shift is inconsistent with the gauge 
condition. It is invariant, however, under an infmitessimal shift together with 
an infmitessimal ‘corrective’ gauge transformation that restores the gauge fix¬ 
ing 


G(x) = ^1 - -r)(x) ■ crj ; U s (x o) = ^1 - ^e 3 (xo)a 3 j , 

Using the invariance of the measure under combination of a shift and a cor¬ 
rective gauge transformation we obtain 


d 


E 


dr)k(z 


d 


Se^Oo) jr? de n( z o) dr) k (z) 


Zw = 0 . 


In shorthand notation 0, the Ehrenfest relation reads 


(W 3 ), 


(W 3 ), 


+ W 3 


(A fp) a 


A FP 


+ 


(A F?)f 


A FP 


+ /3(5)J )=0. (1) 


Gauge fixing has introduced three new terms: 


(W 3 \ 


comes from the corrective gauge transformation acting on the 


source which is U(l) invariant but not SU(2) invariant. 


(A F p) 


(A fp)u 


is the effect of the shift on the Faddeev-Popov determinant, 
is due to the corrective gauge transformation of the Faddeev- 


Popov determinant. 
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The last two derivatives are subtle. The key is to first consider the constraint 
up to first order in the shift and the corrective gauge transformations. Imposing 
that it is still zero fixes the {??}. 


Fi(x) + 


dFj(x) 

de^zo) 


c m( 2 o) + 

k,z 


dFj(x) 

drik(z) 


r]k{z) = 0. 


Then dehne the Faddeev-Popov matrix as a derivative of the corrected con¬ 
straint with respect to a general gauge transformation 


Mix;jy T $ 


9gj(y) 


Ff(x) 


9Ff(x) ^ 

de^{z 0 ) 


*0o) 


9Ff (x) , . 


k,z 


dyk(z 


9=0 


Finally we evaluate the derivative using 

^ =Tr[M-\M),]. 

A check of the Ehrenfest theorem is given in Table 1. Some of the individual 
terms on the right hand side require a 2N x 2 N matrix inversion, where N is 
the lattice volume. Hence to test the result numerically, we chose as small a 
lattice as possible. The exactness of the theorem does not involve the size of 
the lattice which is 4 4 in Table 1. The last column employs a different source. 
The links making up the plaquette are replaced by the diagonal parts only as 
a second test of the theorem. 

Again by separating the links into diagonal and off-diagonal parts we get 
the final form of the Ehrenfest-Maxwell relation. 


= (4 vn -) + Jt tatic 


+ JI 




The right hand side consists of a sum of conserved currents. The first term 
comes from the excitation of the charged coset fields, the static term has an 
extra non-local contribution coming from the corrective gauge transformation, 
and the last two contributions are from the ghost fields. These terms give 
a non vanishing charge density cloud around a static source. The lefthand 
side can be used as a lattice operator to measure the total charge density and 
does not require the matrix inversions needed to measure the individual terms 
separately that limited the numerical tests to small lattices. 
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Source: 

Ehrenfest term 


W 3 (U -> D) 


W 3 


0.65468(10) 0.63069(20) 


0.06095(7) 0.04463(4) 


0.00127(21) 0.00132(50) 


0.00529(3) 0.00564(3) 

-0.72246(68) -0.68275(50) 

-0.00026(77) -0.00045(647 


Tabic 1: Terms in the Ehrenfest relation, Eqn.(l), on a 4 4 lattice at /3 = 2.5. The column 
labeled W 3 corresponds to the source described in the text. In the second column the source 
links are replaced by their diagonal parts of the links to test a second source. The theorem 
gives zero for the sum. 

As a simple application we use this definition of flux to calculate divE on 
a source and the total flux away from the source. In Table 2, we see that the 
total integrated flux on a plane between the charges plus the flux on a back 
plane of the torus is larger than the divE on the source. The interpretation 
is the bare charge is dressed with same polarity charge by the interactions 
and the neighborhood has a cloud of like charge. Hence there is antiscreening. 
This charge density has contributions from all terms in the Ehrenfest relation. 
Table 2 also shows that the interactions increase the charge on the source itself. 

3 Summary 

We have seen that the electric charge induced by Abelian Wilson loop must 
be reinterpreted. The coset fields renormalise the charge of the loop as mea¬ 
sured by |(A“F„ /i }| and charge is also induced in the surrounding vacuum. 



Zero 
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0 

divE(cl.pt. charge) 
i 

~ 3 

divE (on source) 

total flux 

10.0 

0.1 

0.1042(1) 

0.0910(8) (mid) 

(almost 



0.0148(8) (back) 

classical) 



0.1092(8) (total) 

2.4 

0.4166 

0.5385(19) 

0.7455(70) (mid) 
0.0359(72) (back) 
0.7815(95) (total) 


Table 2: divE = (A„ V y \). normalized to 4 for a ‘classical’ point charge, measured on a 3 X 3 
Wilson loop source on an 8 4 lattice. Integrated electric flux is measured on the midplane 
centered on the Wilson loop and on a plane on the far side of the torus, and the sum being 
the total flux. 


Full SU(2) has antiscreening/asymptotic freedom of color charge, and in the 
maximal Abelian gauge alone have we seen analogous behaviour, in that the 
source charge is increased and induces charge of like polarity in the neighboring 
vacuum. Whether this renormalization of charge can account for the reduction 
of the string tension upon Abelian projection in this gauge is not clear. The 
improved field strength expression defined by the Ehrenfest identity does-nrot 
coincide with the lattice versiontdl of the’t Hooft field strength operatorli3. 


Acknowledgments 

We wish to thank G. Di Cecio, A. Di Giacomo, G. Gubarev, M.I. Polikarpov for 
discussions. This work was supported in part by United States Department of 
Energy grant DE-FG05-91 ER 40617. Y.S. was supported in part by a Japan 
Ministry of Education Fellowship. 


References 

1. G. ’t Hooft, Nucl. Phys. B 190 (1981) 455. 

2. A.S. Kronfeld, M.L. Laursen, G. Schierholz, U.J. Wiese, Phys. Lett. B 
198 (1987) 516. 

3. T. Suzuki, I. Yotsuyanagi, Phys. Rev. D 42 (1990) 4257. 


7 








4. G.S. Bali, V. Bornyakov, M. Mtiller-Preufiker, K. Schilling, Phys. Rev. 
D 54 (1996) 2863, available as |hep-lat/960301^ . 

5. S. Mandelstam, Phys. Rept. 23 (1976) 245. 

6. J. Stack, S. Neiman, R. Wensley, Phys. Rev. D 50 (1994) 3399, available 
as hep-lat/9404014 . 

7. M.N. Chernodub, M.I. Polikarpov, A.I. Veselov, Nu cl. Phys. B (Proc. 
Suppl.) 49 (199 6) 307, available as hep-lat/951203C . 

8. A. Di Giacomo, hep-lat/ 9802008| , and references therein. 

9. V. Singh, D.A. Browne, R.W. Haymaker, Nucl. Phys. B (Proc. Suppl.) 
30 (1993) 658, available as hep-lat/ 9302010| ; Phys. Lett. B 306 (1993) 
115, available as hep-lat/9301004 ; Y. Matsubara, S. Ejiri, T. Su zuki 
Nucl. Phys . B (Proc. Suppl.) 34 (1994) 176, available as hep- 


lat/9311061 

Y. Peng, R.W. 

available as 

hep-lat/9503015 


10. G.S. Bali, C. Schlichter, K. Schilling, hep-lat/9802005 . 

11. G. Di Cecio, A. Hart, R.W. Haymaker, LSU Preprint No. 269-1998, 


hep-lat/9807001 


12. M. Zach, M Faber, W. Kain z, P. Skala, Phys. Lett. B 358 (1995) 325, 
available as hep-lat/9508017. 

13. G. Di Cecio, A. Hart, R.W. Haymaker, Nucl. Phys. B (Proc. Suppl.) 
63 A-C (1998) 525, available as |hep-lat/9709084 . 

14. K. Bernstein, G. Di Cecio, R.W. Haymaker, Phys. Rev. D 55 (1997) 
6730, available as |hep-lat/9606018. 

15. G. ’t Hooft, Nucl. Phys. B 79 (1974) 276. 





























